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The  membraneless  microfluidic  fuel  cell  (MFC)  is  a  promising  micro-scale  power  source  with  potentially 
wide  applications.  MFC  commonly  relies  on  the  co-laminar  microfluidic  platform  in  which  redox  streams 
flow  in  parallel  in  a  microchannel.  The  nature  of  this  cell  architecture  limits  the  mass  transport  inside  the 
cell,  often  resulting  in  low  power  density.  To  overcome  the  issues,  we  propose  an  innovative  concept  of 
chaotic  flow-based  fuel  cell  (CFFC),  which  is  built  on  a  counter-flow  microfluidic  platform  with  the  flow 
channel  patterned  with  micro-ridges.  A  CFD/electrochemical  model  is  used  to  predict  the  performance 
and  investigate  the  underlay  mechanism  of  the  CFFC.  Two  theoretical  upper  bounds,  i.e.,  the  limiting 
current  and  limiting  fuel  conversion  for  conventional  MFC,  are  derived.  Through  the  results,  it  is  found 
that  the  generation  of  chaotic  flow  inside  the  patterned  activation  zone  enables  the  CFFC  to  exceed  the 
theoretical  limitations  and  work  with  over-limiting  current  for  high-power  output.  Meanwhile,  the  inter¬ 
facial  mixing  and  crossover  is  minimized  by  the  counter-flow  microfluidics,  allowing  for  over-limiting 
fuel  conversion  to  useful  electricity  output.  The  achievement  of  unprecedented  operating  regime  demon¬ 
strated  in  this  study  open  up  a  new  direction  towards  optimization,  operating  and  design  of  the  MFC. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  membraneless  microfluidic  fuel  cell  (MFC)  has  recently 
drawn  much  research  interest  as  a  promising  micro-scale  power 
source  for  on-chip  devices  and  microelectronics.  The  simple  struc¬ 
ture  and  working  principle  make  MFC  competitive  among  other 
types  of  micro-fuel  cells,  e.g.,  micro  PEMFC  and  DMFC  [1,2].  The 
original  design  of  MFC,  which  is  still  widely  adopted  in  recent  pub¬ 
lished  works  [2-6],  relies  on  the  co-laminar  microfluidic  platform  in 
which  redox  streams  flow  in  parallel  in  a  microchannel  (typically  in 
a  Y-shape)  with  anode  and  cathode  on  the  opposite  channel  walls. 
In  the  laminar  regime,  mass  transport  and  fuel-oxidant  mixing  is 
limited  to  diffusive  effects;  thus  restricted  to  the  interfacial  width. 
Since  the  two  redox  streams  are  naturally  separated,  the  need  for 
a  membrane  is  eliminated.  Meanwhile,  protons  can  still  efficiently 
transport  between  the  two  electrodes  through  electrolyte  (Fig.  la). 
Such  kind  of  MFC  is  commonly  specified  as  laminar  flow-based  fuel 
cell  (LFFC)  [5,6].  Reviews  on  the  working  principles,  designs  and 
recent  development  of  LFFC  can  be  found  elsewhere  [2,4]. 
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To  optimize  the  cell  performance,  one  must  minimize  the  mass 
transfer  at  the  liquid-liquid  interface  to  prevent  the  fuel  crossover, 
while  at  the  same  time,  mass  transfer  through  the  electrode  con¬ 
centration  boundary  layer  should  be  enhanced  for  higher  current 
density  and  fuel  utilization.  These  two  conflicting  requirements 
pose  enormous  challenges  to  the  manipulation  of  the  transport 
patterns  at  such  a  small  scale.  For  example,  simply  increase  in  the 
inlet  flow  rate  may  improve  the  cell  current  density  and  eliminate 
the  fuel  crossover.  However,  the  fuel  utilization  will  decline  due 
to  insufficient  residence  time  for  surface  reactions.  To  date,  seek¬ 
ing  for  an  effective  solution  to  control  and  optimize  the  transport 
behaviors  in  MFC  is  still  an  outstanding  task. 

Shear-driven  chaotic  flows  have  been  noticed  for  their  effects 
on  the  improvement  of  mass  transport  in  microchannels  [7-9]. 
Application  to  the  MFCs  can  be  found  in  Yoon  et  al.  [10],  report¬ 
ing  an  increase  in  the  fuel  conversion  rate  up  to  40%  by  patterning 
micro-ridges  at  the  bottom  of  LFFC  channel.  However,  the  trans¬ 
verse  secondary  flows  generated  by  the  ridges  may,  at  the  same 
time,  intensively  mix  the  two  streams  within  a  short  residence  time 
[7,8].  Consequently,  fatal  LFFC  operation  will  easily  occur  due  to  the 
sharp  increase  in  parasitic  current.  However,  in  the  experiment  by 
Yoon  et  al.  [1 0],  the  fuel  crossover  and  parasitic  effects  were  ignored 
by  using  a  model  redox  of  ferricyanide/ferrocyanide  instead  of  a 
real  fuel/oxidant  couple.  Therefore,  their  results  may  not  be  fully 
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Nomenclature 

Acell  interfacial  area  of  the  computational  cell  (m2 ) 
c  concentration  (M) 

C  fuel  conversion  rate 

Da  Damkohler  number 

Dj  diffusivity  of  species  i  (m2  s-1 ) 

F  Faraday  constant,  96,487  (C  mol-1 ) 

h  channel  height  (m) 

/  current  (A) 

j,  jo  current  density,  exchange  current  density  (A  irr 2 ) 

L  length  of  the  electrode  (m) 

Mi  molar  mass  of  species  i  (kg  mol-1 ) 
n  number  of  electrons  exchanged  at  the  rate-limiting 

step 

nt  number  of  total  electrons  exchanged 

p  pressure  (Pa) 

R  universal  gas  constant  (8.314  J  mol-1  K-1 ) 

Si  source  term  of  specie  i  (kg  s-1  m-3 ) 

source  term  in  charge  conservation  (V  s-1  m-3 ) 

T  temperature  (I<) 

U,  Utrans  velocity,  transverse  slip  velocity  (m  s-1 ) 

Vceu  volume  of  the  computational  cell  (m3 ) 

Vop  operating  voltage  (V) 

w  channel  width  (m) 

Yi  mass  fraction  of  species  i 

Greek  letters 

a  charge-transfer  coefficient 

r\  overpotential  (V) 

p  density  of  the  reactant  streams  (kg  m-3 ) 

0  potential  (V) 

a  conductivity  (S  nrr 1 ) 

r  stress  tensor  (Pa) 

X  reaction  order 

®b  normalized  bulk  concentration 

Subscripts 

a  anode  reaction 

c  cathode  reaction 

/  formic  acid 

o  oxygen 

lim  limiting 

0  standard  or  reference 


applicable  to  a  practical  LFFC.  To  overcome  the  interfacial  mixing 
problem  while  maintaining  enhanced  mass  transport  to  the  elec¬ 
trode,  we  propose  a  novel  concept  of  chaotic-flow  based  fuel  cell 
(CFFC)  that  is  built  on  a  counter-flow  microfluidic  platform.  In  the 
counter-flow  network,  the  activation  zone  (vicinity  of  electrode  and 
micro-ridges)  and  mixing  zone  are  spatio-temporally  separated  to 
facilitate  control  of  the  transport  behaviors  in  microfluidics  in  order 
to  fulfill  different  requirements  on  each  functional  region. 

The  CFFC  is  a  brand  new  research  area  that  is  distinct  from  LFFC 
and  any  other  types  of  fuel  cells,  and  thus,  need  to  be  further  investi¬ 
gated  from  device-level  performance  characteristics  to  microscopic 
level  governing  mechanisms.  Further  development  for  optimiza¬ 
tion  of  the  cell  design  and  operation  is  also  urgent.  Mathematical 
models  have  proven  important  in  understanding  and  optimizing 
a  new  microfluidic  device,  such  as  CFFC.  We  have  successfully 
developed  a  CFD/electrochemical  model  and  used  the  model  to 
investigate  the  micro-scale  physico-electrochemical  interactions  in 
MFCs  [  1 1  -1 3  ].  In  this  study,  we  extend  and  improve  the  model  to  be 
capable  of  modeling  CFFC  of  air-breathing  cathode  design.  The  LFFC 


built  on  counter-flow  microfluidic  network,  as  reported  in  a  recent 
communication  [14],  is  also  included  in  this  study  as  a  contrast 
for  comparison.  Through  the  modeling  results,  we  predict  that  the 
CFFC  can  be  operated  at  an  over-limiting  current  condition,  which  is 
not  possible  for  LFFC.  Unprecedented  current  density  and  fuel  con¬ 
version  rate  are  predicted  to  exceed  the  theoretical  upper  bounds 
for  conventional  MFC.  The  results  suggest  that  CFFC  is  a  promising 
innovation  which  can  contribute  to  the  future  development  of  MFC. 

2.  Methodology 

2.1.  Numerical  model 

2.1  .1.  Physical  problem  and  computational  domain 

An  illustration  of  the  CFFC  considered  in  this  study  is  shown  in 
Fig.  lb.  The  CFFC  consists  of  an  anodic  activation  region,  a  cathodic 
activation  region  and  a  mixing  region.  The  fuel  and  cathodic  sup¬ 
porting  electrolyte  are  fed  through  opposite  inlets  of  the  anode  and 
cathode,  respectively.  The  two  streams  meet  at  the  mixing  region 
and  the  ionic  charge  transport  takes  place  at  the  interface.  Planar 
electrodes  are  assembled  as  the  top  boundaries  of  the  cell.  The  cath¬ 
ode  is  a  gas  diffusion  electrode  that  acquires  oxygen  from  outside 
the  cell.  As  shown  in  Fig.  la,  the  charge  transport  pattern  in  CFFC 
is  different  from  that  in  LLFC,  but  similar  to  that  in  a  MFC  using 
graphite  rods  as  the  electrodes  [15].  To  generate  the  chaotic  flow, 
micro-ridges  are  patterned  at  the  bottom  of  the  cell  to  produce 
rotating  eddies  and  vortices  in  the  flow.  The  ridges  are  designed  to 
repeating  asymmetric  patterns.  As  such,  Lagrangian  chaos  will  be 
generated  [16,17].  Therefore,  such  microfluidic  fuel  cell  is  named 
as  chaotic  flow-based  fuel  cell. 

The  computational  domain  is  built  accordingly.  As  suggested 
by  previous  studies  [8,16],  the  chaotic  effect  of  the  micro-ridges  on 
the  flow  can  be  mimicked  by  replacing  the  detailed  geometry  of  the 
ridges  with  a  flat  boundary  with  transverse  slip  velocity,  as  shown 
in  Fig.  1  c.  The  relationship  between  the  slip  velocity  magnitude  and 
the  geometric  parameters  of  the  ridges  is  well  stated  elsewhere 

[8.18] .  This  simplification  on  the  computational  domain  will  not 
sacrifice  the  accuracy  of  the  transport  and  chemical  behaviors  in 
the  bulk  flow  and  to  the  electrodes,  which  is  the  main  focus  in  the 
study,  although  it  may  simplify  the  details  on  the  patterned  surface 
[16].  The  treatment  of  slip  boundaries  has  been  justified  elsewhere 

[8.16.18] . 

2.1.2.  Electrochemical  model 

Considering  an  air-breathing  CFFC  running  on  formic  acid,  the 
electrochemical  system  becomes 

Anodereaction  :  2HCOOH  2C02  +  4H+  +  4e~ ;  ( 1 ) 

Cathodereaction  :  02  +4FI+  +4e-  ->  2FI2O.  (2) 


For  both  electrodes,  their  reverse  reactions  are  negligible  during 
MFC  operation  [  1 9].  Therefore,  their  reaction  rates  can  be  expressed 
by  Tafel  law. 


In  the  above  equations,  there  are  several  unknown  variables,  i.e., 
fuel  and  oxidant  concentration  at  the  electrode  surface  and  elec¬ 
trode  overpotential,  which  can  be  derived  from  subsequent  CFD 
model.  Particularly,  the  oxygen  concentration  at  the  electrode  sur¬ 
face  is  given  as  a  constant  equivalent  to  its  fraction  in  the  air.  It  is 
because  a  recent  study  [20]  revealed  that  the  oxygen  concentra- 
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Fig.  1.  Illustration  of  CFFC:  (a)  working  principles  and  charge  transport  in  co-laminar  and  counter-flow  MFCs,  (b)  detailed  geometry,  and  (c)  approximation  adopting  flat 
boundary  with  transverse  slip  velocities. 


tion  polarization  at  the  cathode  of  MFC  is  negligible,  providing  the 
external  oxygen  supply  to  the  cathode  is  sufficient. 

2 A. 3.  CFD  model 

In  the  CFD  model,  the  coupled  continuity  equation, 
Navier-Stokes  equations,  species  transport  equation  and  electric 
potential  equation  are  solved  in  the  3D  domain. 


V  .  ( pU )  =  0 

(5) 

V  .  ( pUU )  =  -Vp  +  V  ■  r 

(6) 

V  ■  (pUYi)  -  V  •  (pD,VYf)  =  S,- 

(7) 

aW2<p  =  S# 

(8) 

The  species  transport  equation  is  solved  for  formic  acid.  The 
mass  balance  is  ensured  by  assuming  another  species  of  “bulk 
mixture”,  whose  mass  fraction  is  set  to  be  unity  minus  the  mass 
fractions  of  formic  acid.  The  source/sink  terms,  S,  and  5^,  cor¬ 
respond  to  the  volumetric-basis  species  depletion  and  current 
production  due  to  the  electrochemical  reactions,  respectively.  They 
are  confined  in  a  layer  of  control  volumes  neighboring  the  elec¬ 
trode/electrolyte  interface  and  zero  in  other  computational  zones. 


_  _  MjAce ii 

n^JVceii* 


Sj)  = 


^cell ; 

Vcel* 


O) 

(10) 


2 A  A.  Boundary  conditions 

The  boundary  conditions  of  the  above  equations  include  non¬ 
slip,  non-permeable  conditions  at  channel  walls  and  electrodes, 
designated  slip  velocity  for  mimic  ridged  boundaries,  designated 
velocities  and  concentrations  at  inlets,  and  the  pressure  outlet, 
etc.  Note  that  zero-flux  boundaries  are  specified  to  the  elec¬ 
trodes  here  for  the  species  because  any  production  or  consumption 
has  been  included  as  a  source  term  appearing  in  the  governing 
equations. 


2. 1 .5.  Numerical  procedures 

Input  parameters  of  the  model  are  summarized  in  Table  1, 
which  have  been  adopted  in  our  previous  studies  [13].  Finite 
Volume  Method  code  FLUENT  6.3  [21]  is  used  to  solve  the  equa¬ 
tions  numerically.  To  minimize  the  numerical  diffusion,  structural 
meshes  under  strict  grid  independence  control,  and  higher  order 
discretization  schemes  are  used,  i.e.,  the  semi-implicit  method  for 
pressure  linked  equations  (SIMPLE)  algorithm  with  the  quadratic 
upwind  interpolation  of  convective  kinematics  (QUICK)  scheme, 
adopted  for  the  pressure-velocity  coupling  and  discretization 
schemes. 


2.2.  Limiting  current  and  fuel  conversion 

In  the  analysis,  dimensionless  velocity  and  current  are  intro¬ 
duced,  enabling  the  comparison  of  microfluidic  fuel  cells  with 
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Table  1 

Key  input  parameters  [13]. 


Parameter 

Value 

Electrode  kinetics 

da 

0.5 

Ole 

0.497 

ja,o  (Am-2) 

0.1194 

jc, o  (Am-2) 

3.1  x  10“7 

na 

1 

nc 

1 

Xa 

1 

Xc 

1 

Fluidic  properties 

p  (kg ITT3) 

1008.6 

li  (kgirr1  s-1) 

0.001  (water) 

o  (Sm-1) 

0.4 

Df(  m2  s_1) 

5xlO-10 

Inlet  conditions 

C/(  M) 

1.05 

T(K) 

297 

Fig.  2.  Comparison  between  experimental  data  [3]  and  model  results  on  the  cell 
performance. 


different  operating  conditions  and  geometries.  They  are  defined 

We  further  discuss  the  fuel  conversion  rate  as  a  function  of  U * 

in  a  LFFC.  The  hydrodynamics  in  the  cell  is  pressure-driven  flow. 

Therefore,  the  bulk  concentration  of  the  fuel  when  the  flow  leaves 

,  the  activation  region,  0b,  can  be  analytically  solved  [23]. 

and 


r  =  1 


2  w 


ntFhLcinDi 


U*  =  U 


4  w2 
leT 


(12) 


The  limiting  current  of  a  fuel  cell  is  usually  defined  as  the  point 
in  the  polarization  curve  where  its  shape  becomes  vertical,  and  no 
further  increase  in  current  can  be  obtained  regardless  of  the  applied 
overpotential,  due  to  insufficient  mass  transport  rate  compared  to 
the  electrode  reaction  rate.  In  LFFC,  the  limiting  current  correlates 
with  inlet  velocity  [15,22], 


^lim 


lu*  (U*  <  6.859) 

1.849 U*V3  ([/*>  6.859) 


(13) 


0b  =  exp 


(14) 


where  the  first  eigenvalue,  ,  is  determined  numerically  and  fitted 
approximately  in  the  form  of  [23] 

^  =  1.4304(1  -  exp(-0.68Da))  +  exp  (15) 

At  the  critical  condition  when  /  =  /lim,  the  Damkohler  number 
(Da)  for  anode  becomes  infinity.  At  that  condition,  <9 b  reaches  its 


Fig.  3.  Comparison  of  fuel  distribution  (left,  unit:  M)  and  flow  patterns  (right,  colored  by  the  coordinates  of  the  origins  of  the  patterns)  for  CFFCs  based  on  (a)  counter-flow 
platform  and  (b)  co-laminar  flow  platform. 


J.  Xuan  et  al.  /  Journal  of  Power  Sources  196(2011)  9391-9397 


9395 


Fig.  4.  Comparison  of  flow  patterns  (left,  colored  by  fuel  concentration  at  local  position)  and  fuel  distribution  (right,  unit:  M)  for  CFFC  with  Utrans/U  equal  to  (a)  0,  (b)  0.5,  and 
(c)  1. 


minimum  value  and  the  corresponding  limiting  fuel  conversion 
rate,  Qim,  can  be  expressed  in  the  way  of 

Ciim  =  1  —  =  1  -  exp  (— ^— )  (16) 

So  far,  we  have  reviewed  and  derived  two  theoretical  upper 
bounds  for  the  operation  of  LFFC,  i.e.,  Jlim  and  Clim.  Based  on 
that,  two  operating  conditions  of  a  microfluidic  fuel  cell  can  be 
defined,  which  are  under-limiting  current  condition  (7</lim  and 
C<Clim)  and  over-limiting  current  condition  (/>/jim  and  C>Clim). 
Normally,  for  practical  operation  of  LFFC,  various  losses  lead  to  the 
cell  performance  deviating  from  the  theoretical  upper  bounds  to  an 
under-limiting  current  condition.  In  the  subsequent  sections,  we 
demonstrate  that  the  operation  of  CFFC  outperforms  the  theoreti¬ 
cal  upper  bounds  posed  on  LFFC  and  unprecedented  high  current 
density  and  fuel  conversion  are  obtained. 

3.  Results  and  discussion 

3.1.  Model  validation 

Validation  of  the  model  is  conducted  by  comparing  the  indi¬ 
vidual  anode  and  cathode  performances  of  an  air-breathing  LFFC 
obtained  from  the  model  and  experimental  results  [3].  The  geo¬ 


metric  and  operating  input  parameters  of  the  simulation  follow 
the  experimental  study  by  Jayashree  et  al.  [3].  Fig.  2  shows  the 
results  of  the  comparison,  from  which,  reasonable  agreement  can 
be  found.  Moreover,  implementations  of  the  validation  of  the 
present  model  have  also  been  conducted  several  times  in  our  pre¬ 
vious  works.  We  have  shown  that  the  model  is  able  to  accurately 
predict  the  convection-diffusion-reaction  interactions  [11,12]  and 
electrochemical  performance  [12,13]  of  microfluidic  fuel  cells. 

3.2.  Interfacial  transport 

As  mentioned  earlier,  the  chaotic  flows  in  the  MFC  may  mix 
the  anodic  and  cathodic  streams  and  cause  fuel  crossover  problem. 
Therefore,  it  is  important  to  evaluate  the  interfacial  mixing  patterns 
of  the  CFFC  to  prove  the  feasibility  of  the  concept.  Fig.  3  compares 
the  hydrodynamic  conditions  and  interfacial  transport  patterns  in 
CFFCs  based  on  counter-flow  and  co-laminar  flow  platforms  under 
a  velocity  Utrans/Ll=0.5.  For  clear  illustration  and  explanation,  the 
electrode  reactions  are  not  included  at  the  current  stage  (i.e.,  7  =  0). 
It  is  found  that  for  CFFC  based  on  counter-flow  network,  the  inter¬ 
facial  mixing  is  well-controlled.  Within  all  cases  studied  in  this 
paper,  no  crossover  of  the  fuel  is  observed.  It  is  because  the  trans¬ 
verse  secondary  flows  generated  by  micro-ridges  are  parallel  to 
the  interface  of  the  two  streams,  and,  thus,  will  not  contribute  to 
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Utrans/U 


0.1  1  10  100  1000 


u* 


Fig.  6.  (a)  Dimensionless  current  and  (b)  fuel  conversion  versus  U. 


the  interfacial  mixing.  Moreover,  due  to  the  dominating  viscous 
effects,  the  chaotic  flow  will  become  a  steady  laminar  flow  at  the 
mixing  region,  once  it  leaves  the  patterned  activation  region.  This 
hydrodynamic  behavior  benefits  the  minimization  of  the  interfacial 
mixing  by  eliminating  the  harmful  convective  flux  at  the  stream 
interface.  On  the  contrary,  CFFC  based  on  co-laminar  flow  network 
shows  serious  mixing  and  fuel  crossover  problems  inside  the  chan¬ 
nel,  which  is  similar  with  a  microfluidic  mixer  [16,17].  The  results 
indicate  that  the  counter-flow  network  is  potentially  a  much  more 
promising  platform  to  realize  the  CFFC  concept. 

3.3.  Over-limiting  current  behavior 

We  consider  a  CFFC  with  1  =  5  mm,  w  =  0.5  mm,  U  =  1  mm  s-1 ,  and 
Vop  =  0.1  V.  For  air-breathing  MFC,  the  cell  current  density  and  fuel 
conversion  is  primarily  determined  by  the  anode  part  [5,10,24]. 
Fig.  4  shows  the  flow  patterns  and  fuel  concentration  field  in 
the  anode  activation  region  of  the  CFFC  at  Utrans/O  =  0,  0.5  and  1, 
referring  to  laminar,  low  chaotic  and  high  chaotic  flow  regimes, 
respectively.  It  is  found  that  at  dtrans/^ = 0  (i.e.,  equivalent  to  a  LFFC), 
the  flow  in  the  channel  is  laminar  and  the  cell  is  working  at  lim¬ 
iting  current  condition  with  a  gradually  developed  concentration 
boundary  layer  adjunct  to  the  electrode.  With  the  increase  in  the 
FftransM  the  secondary  helical  flow  will  be  generated  by  the  pat¬ 
terned  ridges  and  the  flow  in  the  activation  zone  becomes  more 
chaotic  (Fig.  4  left  column).  At  the  same  time,  mass  transfer  to  the 
electrode  is  enhanced  in  the  channel,  leading  to  steeper  concentra¬ 
tion  variation  in  the  boundary  layer  and  low  fraction  of  the  fuel  in 


the  bulk  (Fig.  4  right  column).  Therefore,  as  shown  in  Fig.  5,  I*  and 
C  are  enhanced  with  the  increase  in  the  L/trans/ff,  resulting  in  an 
over-limiting  current  condition  in  which  />/lim  and  C>  Clim.  Essen¬ 
tially,  the  performance  exceeding  theoretical  upper  bounds  can  be 
attributed  to  the  benefit  of  the  chaotic  flow  regime  in  the  CFFC, 
since  the  two  theoretical  upper  bounds  are  calculated  based  on 
pressure-driven  laminar  flow  condition.  Note  that  the  current  den¬ 
sity  and  fuel  conversion  rate  can  be  simultaneously  enhanced  by 
increasing  f/trans/ff  by  proper  design  of  the  micro-ridge  patterning, 
which  does  not  suffer  from  the  common  trade-off  limitation  in  LFFC 
optimization  [24]. 

We  then  examine  the  effect  of  operating  flow  rate  on  the  per¬ 
formance  of  the  CFFC.  In  the  model,  the  design  and  operating 
parameters  of  the  CFFC  are  set  to  be  identical  as  the  previous 
case,  except  that  Utrans/Ff  is  fixed  to  1  and  U*  becomes  a  variable. 
A  LFFC  based  on  the  counter-flow  microfluidic  platform  with  the 
same  geometry  and  operation  condition  but  without  patterning 
is  included  for  comparison.  Fig.  6a  and  b  shows  the  change  of  I* 
and  C  against  various  U*t  respectively.  The  over-limiting  current 
and  under-limiting  current  regions  are  distinguished  by  identify¬ 
ing  J*lim  and  Clim  according  to  Eqs.  (13)  and  (16).  It  is  observed  that, 
generally,  the  current  density  increases  while  the  fuel  conversion 
decreases  with  increasing  inlet  velocity  for  both  the  CFFC  and  LFFC. 
Flowever,  the  CFFC  could  exceed  the  critical  values  of  J*lim  and  Clim 
to  work  at  over-limiting  current  region  in  a  wide  range  of  operating 
flow  rate,  while  the  LFFC  can  only  be  operated  at  under-limiting  cur¬ 
rent  region.  The  advantage  of  CFFC  is  clearly  shown.  At  higher  flow 
rate  condition,  both  curves  for  CFFC  and  LFFC  return  to  the  under¬ 
limiting  current  region.  It  is  because  the  transport  limitation  is  not 
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reached  with  high  inlet  velocity:  high  flow  rate  leads  to  high  mass 
transport  to  the  electrode;  when  the  mass  transport  rate  becomes 
faster  than  the  surface  reaction  rate,  electrochemical  and  electri¬ 
cal  aspects  of  the  cell  will  play  the  primary  role  to  determine  the 
cell  performance,  i.e.,  non-ideal  reaction  kinetics  and  ohmic  resis¬ 
tance  becomes  significant  [15].  It  indicates  that  using  chaotic  flow 
to  enhance  the  performance  of  MFC  will  show  distinct  effect  when 
the  mass  transport  in  the  cell  is  critical,  i.e.,  at  high  current  con¬ 
dition  and  low  flow  rate.  On  the  contrary,  at  activation  and  ohmic 
regime  of  cell  operation,  the  strategy  becomes  less  significant.  It 
practice,  to  approach  ultra  high  fuel  utilization  (i.e.,  near  100%),  the 
cell  must  be  operated  at  low  flow  rate  and  high  current  condition, 
under  which  condition  the  CFFC  technique  would  be  particularly 
useful. 

Finally,  the  J-V  characteristics  of  the  CFFC  and  LFFC  on  counter¬ 
flow  network  are  predicted  and  compared  in  Fig.  7.  It  is  found  that 
at  high  current  operating  condition,  the  performance  of  LFFC  is  lim¬ 
ited  by  mass  transport  to  the  anode,  resulting  in  limiting  current 
behavior.  It  means  that  the  cell  current  is  restrained  by  the  theo¬ 
retical  upper  bound,  Jlim,  and  becomes  irrelevant  to  the  operating 
potential.  Similarly,  the  J-V  curve  for  CFFC  yields  the  same  trend 
as  the  one  for  LFFC  at  low  to  medium  current  conditions.  Yet,  the 
improved  performance  of  CFFC  becomes  obvious  at  high  current 
condition  when  the  concentration  overpotential  turns  significant. 
Furthermore,  the  operating  condition  extends  to  the  over-limiting 
current  region  without  the  restrain  by  the  theoretical  limitation, 
i.e.,  /lim.  From  the  results,  it  is  calculated  that  the  CFFC  design 
enables  increase  of  maximum  current  density  from  137  mA  cm-2 
to  181  mAcrrr2,  and  maximum  power  density  from  30  mW  cm  2 
to  35  mW  cm  2. 

However,  both  cell  potentials  of  LFFC  and  CFFC  drop  rapidly  with 
the  increase  in  the  current  density  at  low  to  medium  current  oper¬ 
ating  condition,  which  can  be  mainly  attributed  to  the  high  ohmic 
overpotential  involved.  It  is  because  the  way  of  arrangement  of  two 
electrodes  in  counter-flow  platform  causes  longer  route  for  the 
charged  ions  to  transfer  between  two  electrodes  compared  with 
its  co-laminar  flow  counterpart  (Fig.  la).  To  minimize  the  ohmic 
loss,  the  electrode  length  shall  be  reduced  and  cell  scale  shall  be 
miniaturized  accordingly.  A  decrease  in  the  electrode  length  in  the 
cell  design  can  lead  to  an  increase  in  the  cell  current  density  since 
most  of  the  current  are  collected  within  the  region  near  the  inlet 
[24].  However,  the  fuel  utilization  will  be  lowered  due  to  the  reduc¬ 
tion  in  the  residence  time  of  reaction.  Since  the  fuel  conversion 
rate  can  be  compensated  by  adopting  higher  Utrans/U  or  lower  U 
in  CFFC  operation,  the  issue  of  low  fuel  utilization  would  be  partly 
overcome. 


4.  Conclusion 

Manipulation  of  the  transport  patterns  in  membraneless 
microfluidic  fuel  cell  is  of  great  importance  to  obtain  high  per¬ 
formance.  Yet,  rare  successful  attempts  have  been  reported.  This 
study  proposed  a  novel  concept  of  chaotic  flow-based  fuel  cell 
built  on  counter-flow  microfluidic  platform  to  fulfill  the  research 
gap.  The  design  benefits  from  the  separation  of  the  activation  zone 
and  mixing  zone,  which  facilitates  to  accurately  manipulate  the 
transport  behaviors  for  each  important  physico-electrochemical 
process.  For  example,  the  mass  transport  at  activation  zone  and 
mixing  zone  will  not  influence  each  other.  A  CFD/electrochemical 
model  is  developed  to  study  the  performance  and  underlay  mech¬ 
anism  of  the  CFFC.  Moreover,  two  theoretical  upper  bounds  (i.e., 
the  limiting  current  and  limiting  fuel  conversion)  for  conven¬ 
tional  MFC  have  been  derived.  We  have  predicted  that  the  CFFC 
can  be  operated  at  the  over-limiting  current  condition  to  exceed 
the  theoretical  upper  bounds  and  achieve  unprecedented  per¬ 
formance  of  high  current  density  and  high  fuel  conversion  rate. 
The  unique  operating  regime  enables  simultaneous  optimiza¬ 
tion  of  current  density,  fuel  utilization  and  interfacial  mixing  by 
proper  design  of  the  micro-ridge  patterning  without  the  trade-off 
limitation.  Similar  achievement  has  never  been  reported  in  the  lit¬ 
erature.  Therefore,  the  study  opens  up  a  new  direction  for  MFC 
research. 
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